Systematic studies of hadron production spectra in collider experiments 
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A shape of invariant differential cross section for hadron production as function of transverse 
momentum is analysed. The systematic analysis of the available data demonstrates a need for a 
modification of the parameterization traditionally used to approximate the measured spectra. The 
properties of the new proposed parameterization are discussed. 



I. INTRODUCTION 



energy (y/s) ranging from 23 to 2360 GeV. 
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There exists a large body of experimental data on 
hadron production in high energy proton-proton, photon- 
proton, photon-photon and heavy ion collisions [TMl5|. 
The spectra of hadrons produced in these collisions are 
characterized by an exponential behavior as function of 
transverse energy (Et) for the bulk of produced hadrons, 
which populate the low Et part of the spectra. This be- 
havior resembles the Boltzmann-like spectrum in classical 
thermodynamics. The exponential shape of the spectra 
changes to a power law for high Et hadrons. This change 
is traditionally interpreted as an onset of the perturba- 
tive QCD regime of hadron production. These features of 
the spectra shape are found to be universal for any type 
of colliding particles. Therefore, it is tempting to find 
one universal smooth functional form, which describes 
the spectra of produced hadrons in the whole available 
Et range for different energies and types of colliding par- 
ticles. The parameters of such universal functional form 
are expected to vary for different collision energies, types 
of colliding particles and types of produced hadrons. A 
study of the variations of these parameters provides a 
unique information on the hadron production dynamics. 

In the present paper the experimentally measured in- 
clusive spectra of long-lived charged particles (mainly 
charged pions) produced at central rapidities in the col- 
liding particles center of mass system are analyzed. Fur- 
ther on it is assumed all charged particles being pions 
for simplicity. The present analysis is based on the 
published hadron production measurements made with 
pp— collisions at ISR pQ and LHC [15] . pp— collisions at 
SppS and Tevatron [5J [B], AuAu— collisions with 

different centralities at RHIC [13] , real 77— collisions 
at LEP E] and jp— interactions with different val- 
ues of photon virtuality (Q 2 ) at HERA. The HERA 
data additionally allow to consider two different regimes 
of 7p— interactions: photoproduction at low values of 
Q 2 [3IH1IID] and Deep Inelastic Scattering (DIS) for high 
Q 2 values O [13] . The data for all these inclusive differ- 
ential cross section measurements have been taken with 
a minimum bias trigger conditions and at center of mass 



II. SPECTRUM ANALYSIS 

A typical charged particle spectrum as function of 
transverse energy is shown in Fig [I] 
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FIG. 1. A typical charged particle spectrum fitted using the 
Tsallis-type function Q. 

This spectrum is fitted using the Tsallis-type func- 
tion [To] 
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where Pt is transverse momentum of the produced par- 
ticle, Et — y/ Pj, + M 2 with M equal to the pion mass. 
The parameterization ([!]) has only three free parameters: 
A, T and n. While A is an overall normalization, the T 
and n carry important information on the hadron pro- 
duction dynamics. Since for low Et values the param- 
eterization is reduced to the Boltzmann exponent 
~ exp(—ET/T), the parameter T is a QCD analogy to 
a temperature in classical thermodynamics. Note, this 
analogy isn't straight forward, however. In perturbative 
QCD calculations the value of the parameter n depends 
on the shape of structure functions and partonic con- 
tent of the colliding particles as well as on the differen- 
tial cross sections of parton-parton interactions. As it is 
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FIG. 2. A 
630 GeV g| 



ratio of the data to the Tsallis-type |T])(a,b) and the modified |2|(c,d) fit functions in pp collisions at yfs 
(a,c) and RHIC Au - Au collisions at )/s = 200 GeV/N [H] (b,d). 



demonstrated in Fig [I] the economic Tsallis-type param- 
eterization provides a good overall description of the 
spectrum, thus making the Tsallis-type function broadly 
used to ht the recent measurements at RHIC and LHC. 



III. MODIFICATION OF THE TSALLIS 
FUNCTION 

A. Systematic defects of the fit 

However, a closer look at the fit shown in Fig 1 dis- 
closes systematic defects of the fit. In order to see these 
defects it is convenient to plot a ratio of the data to the 
fit function. Two such ratios are shown in Fig[2^a,b) for 
the data sets provided by the UA1 experiment with pp 
collisions at = 630 GeV g] (Fig and by RHIC 
Au - Au collisions at yTs = 200 GeV/N [14] (Fig $p). 
On both plots one observes dips and bumps around or 
above 1 GeV scale. 

The observed difference between the shapes of data 
spectra and function 0) is also typical for other available 
data sets, not shown here. This observation indicates 



that the true shape of the measured spectra doesn't fol- 
low exactly the Tsallis-type parameterization 0. More- 
over, a larger mismatch between the shapes of the data 
spectrum and fit function is found for J/^f produc- 
tion at the Tevatron |17j . In the following an attempt is 
made to find a function which approximates the shape of 
the inclusive hadron production spectra better than the 
Tsallis parameterization. 



B. Search for the new fit function 

To construct a new fit function a combination of two 
functional forms has been used: the exponential and 
power-law. The both functional forms must be the func- 
tions of scalars Pj, or E^ n = E T — M. This choice 
of these variables is very convenient since the both Pj; 
and E!p n vary from zero to the corresponding kinemati- 
cal limits. Non-true scalars like Pt or Pj. have not been 
considered here. Finally, the parameterization 

= A e cxp (-E^/T e ) + ^— , (2) 

P T dP T (1 + 5^)™ 
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turned out to be in a good agreement with the data. 
The new ratios between data and fit function ^ are 
shown in Fig [5jc,d) for the considered above UA1 and 
RHIC data sets correspondingly. Fig[2jc,d) demonstrates 
a significant improvement of the quality of the spectra 
shape approximation by functional) with respect to that 
using the Tsallis-type function (llj). In the following we 
present the arguments why this improvement is not a 
trivial result of increasing the number of free parameters 
in the fit function. 



IV. CORRELATION BETWEEN THE 
PARAMETERS 

The most surprising feature of the new parameteriza- 
tion pi is a strong correlation between the parameters 
T e and T. The dependence of fitted values of the param- 
eter T 2 versus T 2 for charged particles produced in pp, 
pp and Au — Au collisions is shown in Fig [3] 



distribution. The rest bulk of irradiated hadrons repre- 
sents a quasi-thermolized hadronic gas produced with a 
characteristic temperature T e . The spectrum of hadrons 
in this gas has the Boltzman exponential shape. 

A. Contributions of the exponential 
and power law terms 

The contributions of the exponential and power law 
terms of the parameterization ^ to the typical spectrum 
of charge particle produced in pp collisions are shown 
separately on Fig |4j The relative contribution of these 
terms is characterized by a ratio R of the exponential to 
power law terms integrated over P 2 : 

_ A e (2m + 2T e )(n-l) 
U ~ l8AlA-n-T e [6) 

For the spectrum in Fig [4] the exponential term domi- 
nates, while the power law term contributes at the level 
of about 20%. 
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FIG. 3. A correlation between the parameters T 2 and T 2 

This dependence is approximated by a linear function. 
This provides an additional constraint for the parame- 
terization ([2| and therefore reduces a number of free pa- 
rameters used in the fit to the data. Though the physical 
origin of the observed correlation is not quite clear, the 
new constraint helps to minimize uncertainties of the pa- 
rameter values obtained from the fits. This constraint 
will be used further on through this paper. 



V. A TOY MODEL INTERPRETATION 

The form of the parameterization |2]) has a simple toy- 
model interpretation. Within this toy-model a small frac- 
tion of hadrons are produced directly in parton-parton 
interactions of the colliding particles. As required by the 
perturbative QCD the spectrum of particles produced in 
parton-parton interactions is described by a power law 
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FIG. 4. The contributions of the exponential (red) and the 
power law (green) terms to the spectrum of charged particles 
produced in pp collisions. 

The ratio R for the inclusive charged particle spectra 
for pp and pp collisions as function of y/s is shown in 
Fig [5^. Interestingly, this ratio is almost independent of 
the collision energy and equals to about 4. In FigJHJ) the 
ratio R is shown for Au — Au interactions at RHIC as 
function of centrality of the heavy ions collision. As seen 
on Fig [SJd the relative contribution of the exponent term 
reaches minimum values at medium centralities of heavy 
ion collisions. 

On the other hand, the QCD partonic interaction must 
describe any hard scattering process like the heavy quark 
production, or high Q 2 DIS. Indeed, as shown in Fig[6]thc 
spectrum of heavy J /\V quarkonium produced in high en- 
ergy pp collisions at Tevatron [17] follow the pure power 
law distribution and has no room for an exponential term 
of the parameterization pi). In addition, in the high en- 
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FIG. 5. The 
function of yfs 



ratio R of the exponential to power law contributions to the parameterization |2|for pp and pp collisions as 
(a) and for Au — Au interactions at RHIC as function of centrality (b). 




FIG. 6. The spectrum of J/$ produced in high energy pp 
collisions at Tevatron [17] . 
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FIG. 7. The ratio R of the exponential to power law contri- 
butions to the parameterization Q in the high energy DIS, 
photoproduction and 77 collisions. 



ergy DIS, photoroduction and 77 collisions the power law 
term of the new proposed parameterization ([2| dominates 
the produced particle spectra as shown in Fig [7| Thus, 
only the inclusive spectra of charged particles produced 
in pure baryonic collisions require a substantial contri- 
bution of the Boltzman-like exponential term. We found 
it to be particularly interesting that the interactions of 
real photons and protons with real photons (high energy 
photoproduction) have practically no Boltzman-like ther- 
molizcd hadronic final states. 



VI. MAP OF PARAMETERS 

Finally, a map of the parameters T and n for proton- 
(anti)proton, heavy ion, gamma-proton and gamma- 
gamma collision at different energies is displayed in Fig[8] 
There are two clearly distinct trends seen in Fig p] 
The pp and pp collision data show an increase of the 



T— parameter and decrease of the n— parameter with col- 
lision energy y/s increasing. The second trend, where 
the values of both parameters the T and n increase, is 
defined mainly by the RHIC Au — Au collision data at 
yfs = 200 GeV per nucleon. In this case a simultaneous 
increase of the T and n values corresponds to an increase 
of the centrality (or charged multiplicity) of heavy ion 
collisions. Surprisingly, the both trends cross each other 
at medium centralities corresponding to the minimum 
bias Au — Au collisions and pp interactions with energy 
of y/s — 200 GeV. Naively one could expect the single pp 
interaction has more similarity to the very peripheral sin- 
gle nucleon-nucleon interactions. Contrary to that naive 
expectation, the Deep Inelastic Scattering (DIS), jp and 
77 interaction data with yfs ranging approximately from 
100 GeV to 200 GeV belong to the second trend shown 
in Fig [8] and are located on the parameter map (Fig [8]) 
nearby very peripheral heavy ion interactions at about 
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FIG. 8. Map of the parameters T and n for proton-(anti)proton, heavy ion, gamma-proton and gamma-gamma collision at 
different energies. 



the same collision energy per nucleon. 



VII. SHAPE OF THE SPECTRUM AT HIGH P t 



The proposed new parameterization ^ describes the 
data well at low and medium values of the transverse 
momentum of the produced particles. However, as it was 
recently shown by the CDF measurements, in pp inter- 
action the inclusive particle spectrum has a significant 
excess over the simple power law shape for particles pro- 
duced with P T > 10 GeV/c [6]. In Fig [9] it is demon- 
strated, that observed excess is well described by adding 
a second power law function with somewhat lower value 
of the exponent n. This phenomenon, which has no con- 
sistent explanation by now, seems to be found earlier in 



7P and 77 interactions (shown in Fig 10 ) 



Though, in these types of collisions with at least one 
photon involved the excess of the data spectra over a 
single power law function is visible already for particles 
with P T > 3GeV/c. 




[GeV] 

FIG. 9. CDF data [6] fitted by two power law functions with 
different free parameters. 



VIII. CONCLUSION 

In conclusion we have proposed a new parameteriza- 
tion of the spectrum shape of inclusive charged parti- 
cles produced in high energy collisions. This new pa- 
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ramcterization describes the available experimental data 
significantly better than the broadly used Tsallis-type 
parameterization. The proposed parameterization is a 
sum of an exponential (Boltzman-like) and a power law 
(Tsallis-like) terms. The parameters of these two terms 
turned out to be strongly correlated. We observe, that 
the shapes of the power law terms in minimum bias heavy 
ion collisions and in proton- (anti)proton interactions at 
the same collision energy are practically the same. Addi- 
tionally, the shapes of the power law terms in the spectra 
measured with very peripheral heavy ion collisions and 
interactions of high energy particles with photons show 
very close similarity. The difference in size of the ex- 



ponential Boltzman-like contribution to these spectra is 
mainly responsible for the difference in the spectra shapes 
observed in these experiments. 
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